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ABSTRACT 

We report on observations of four southern Anomalous X-ray Pulsars, 
(1RXSJ170849. 0-400910, IE 1048.1-5937, IE 1841-045 and AX J1845-0258), 
obtained at 1.4 GHz using the Parkes radio telescope. Radio pulsations from these 
sources have been searched (i) by directly folding the time series at a number of trial 
periods centered on the value of the spin rate obtained from the X-ray observations; 
(ii) by performing a blind search; (iii) using a code sensitive to single dedispersed 
pulses, in the aim to detect signals similar to those of the recently discovered Rotating 
RAdio Transients. No evidence for radio pulsations with an upper limit of ^ 0.1 mJy 
for any of the four targets has been found. The blind search led to the serendipitous 
discovery of two new pulsars, rotating with a spin period of about 0.7 s and of 92 ms 
respectively, and to the further detection of 18 known pulsars, two of which were also 
detected in the single-pulse search. 

Key v^rords: pulsars: searches - AXPs: individual (IE 1048.1-5937, IRXS J170849.0- 
400910, IE 1841-045, AX J1845-0258) 



1 INTRODUCTION 

In the last two decades, a significant amount of ob- 
servational and theoretical effort has been dedicated to- 
wards the understanding of an "unusual" class of pulsars, 
namely the Anomalous X-Ray Pulsars (AXPs) and the Soft 
Gamma-Ray Repeaters (SGRs). These relatively bright X- 
ray sources (lO'^* — 10'^'' ergs~^), discovered either as pecu- 
liar short X-ray bursts or steady X-ray emitters, were soon 
recognised as a distinct class of neutron stars with respect 
to the well known radi o pulsar or X-ray binary populations 
(see [woods & Thompson| [2004l for a review). 

The arguments that lead to this conclusion are the fol- 
lowing: first, AXPs and SGRs are X-ray pulsators with spin 
periods in the range 5-12 s. Furthermore, their rotational en- 
ergy loss E inferred from their spin-down rate [E = I^Q,, 
where I, Q. and Q. are respectively the neutron star mo- 
ment of inertia, angular velocity and angular velocity time 
derivative) is insufficient to power the observed X-ray lu- 
minosity (note that for the large majority of isolated pul- 
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sars, detected both in radio and X-ray band, the X-ray lu- 
minosity is typically 0.1% of the i?: lBecker fc TriimperllQQ'^ : 
IPossenti et al.l2002l) . Second, there is no evidence for a com- 
panion star (through the search for Doppler orbital modu- 
lation of the X-ray s pin period, and through observations 
at other wavelength s; Mcrcg hetti et ai]ll998l : IWilson et alJ 
ll999l:lHulleman et al 200Q.) . which might power their emis- 
sion through the accretion mechanism, at least with a mass 
greater than 0.1 M0. An alternative source of energy is 
therefore required to power the X-ray emission of AXPs and 
SGRs. 

Assuming these sources being isolated X-ray pulsars 
with a dipole magnetic field configuration and purely mag- 
netic dipole losses, their inferred magnetic fields are ex- 
tremely high, B'^ 3.2 X 



' lO''^ G: these 

magnetic fields have been identified as the possible addi- 
tional energy reservoir. In fact, if AXPs and SGRs are 
indeed "magnetars", i. e. isolated highly magnetic neu- 
tron stars, many of their bizarre X-ray properties might 
be explained: the X-ray luminosity is naturally provided 
by the decay of these ultra-strong magnetic fields, while 
their peculiar short bursts may result from magnetic induced 
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cracks on the neutron star surfa ce JPuncan fc ThompsorJ 
ll992l:lThompson fc Duncanll993l: 1995). Moreover, the long 
lasting fl ux enhancements observed in some magnetars 
iRea et al.,..2005 : Mcrcghetti et al, 2005^1 may be associated 
to single fractures on a larger spatial scale, or to the episodic 
onset of a wind energised by frequent, sm all scale fractures 
and / or quasi-steady seismic vibrations iThompson et all 
120021) . In this scenario, the magnetars' magnetic fields are 
all above the so-called "quantum critical field" Be = 4.4 x 
10^"^ G, at which the radio pulsations are expected to be 
suppressed by processes such as photon splitting whic h may 
inhibit pair-production cascades (firing & Harding ll99St l 
responsible for creating the radio emission of the normal ra- 
dio pulsars. However, the discovery of high magnetic fiel d ra- 
dio pulsars (|Camilo et alJ200nl:lMcLaughlin et all2003h and 
transient radio pulsed emission from one AXP llCamilo et alJ 
[2OO6) put many uncertainties with this respect making again 
theoretically plausible a link between AXPs and radio pul- 
sars. 

In the past years, radio observations of AXPs and SGRs, 
performed using different antennas and arrays, have given 
always negative results (e.g. Crawf ord et al^ 200 2) both to 
timing and interferometry searches, but rev ealing in a few 
cases some reliable associations with SNRs (jOaensler et al.l 
I2OOII) and in one case an interesti ng stellar wind bub - 
ble blown by the massive progenitor iG aensler erai]l2005l) . 
However, in some peculiar circumstances radio counter- 
parts were detected: a) the SGRs' Giant Flares, rare, catas- 
trophic and highly energetic events, have been discovered 
powering a radio outburst, which slowly faded with a 



timescale of weeks Frail et al]ll999l : ICameron et al J 1200,4 
iGaensler et ahllioOS l. b) during the outburst of the only 
confirmed transient AXP, XTE J1810-197, the Very Large 
Array (VLA) observations (with angular resolution of 6") re- 
vealed a point source of flux density 4.5±0.5 mjy at 1.4 GHz 
iHalpern et al.ll2005^ . which has recently be en discovered 
as pu lsated at the neutron star spin period dCamilo et alJ 

With this picture in mind we have undertaken a sys- 
tematic deep search for radio pulsations in three con- 
firmed (IRXS J170849.0-400910 , IE 1048.1-5937 , IE 1841- 
045 ) and one candidate AXP (AX J1845-0258 ), visible from 
the southern hemisphere using the Parkes radio telescope. 
Observations and data analysis are described in f|H while 
in il8.1l we report on the obtained results, discussed in the 
context of the magnetar scenario in 



2 OBSERVATIONS AND DATA ANALYSIS 

In Table Q the sample of observed AXPs is presented. For 
each object the name, period and period derivative of X-ray 
pulsations, the derived surface magnetic field, the estimated 
distance D and the association with a SNR or with a nebula 
are listed. 

The radio observations have been performed between 
October 12 an d 15 199 9 using the 13 beams of the multi- 
beam receiver dStavelev- S mith et aLi.l9 96l of the Parkes ra- 
dio telescope (NSW Australia) at a frequency of 1374 MHz, 
with the central beam pointed on the target AXP. In order 
to mitigate the effects of the dispersion of the signals in the 
interstellar medium (ISM) , the total 288 MHz bandwidth is 
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Figure 1. P — P diagram for radio pulsars (dots) and magnetars 
(stars). The dashed line represents the limit for radio emission 
quenching. The four sources surrounded by a square emit in radio 
besides being above this line. The three confirmed AXPs of our 
sample for which a determination of P is available are surrounded 
by a circle. 



splitted into 96 frequency channels each 3-MHz wide. The 
outputs from each channel are summed in polarisation pairs, 
high-pass filtered and 1-bit sampled every 1.0 or 1.2 ms. 
Each source has been pointed for 2.8 hours and for all but 
1E1841— 0450 the observation has been performed twice. 

The data for the central beams have been first analysed 
using the pdm programme: the code takes as input a period 
P and a dispersion measure DM and folds the time series 
according to a number of trial values around the input ones, 
searching for the combination of P and DM for which the 
signal-to- noise ratio S/N is maximised. The period range 
searched for each source has been obtained from the X-ray 
ephemeris and their errors. Since at the time of the data 
analysis most of the X-ray timing solutions had big uncer- 
tainties, the range explored (see Table|5J is much wider than 
would have been necessary with the present knowledge. A 
number of trial DM values, ranging from up to the value 
DMmaa; giviug a maximum pulse broadening in a 3 MHz 
frequency channel of 10% of the pulse period (see Table |2l 
has been explored. HMmax is well beyond what expected 
for a source in the Milky Way on the basis of the available 
models for the distribution of free electrons in the Galaxy 
iTavlor fc CordeJl993l:fCordes fc Lazij2002l) . but may take 
into account the possibility of the presence of dense local 
matter. 

Making use of the relation derived byHhaFeFai] (Eooi) 
between the broadening of the pulse due to the scattering 
caused by the ISM, Stscatt and the DM, we note that already 
a DM ~ 1200 pc cm~'^ results in a broa dening of ~ 10 sec- 
onds at 1.4 GHz. However the relation bv lBhat et al.l (|20^ 
has a very high scattering in principle allowing for reason- 
able values of Stscatt also for much higher DMs. Hence, for 
the relatively fast search performed with pdm we have ex- 
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Name 


P 


P 


B 


D 




association 




(s) 


(10-11 ss-l) 


(1014 G) 


(kpc) 


(pc cm— ^) 




IE 1048.1-5937 


6.456109(5) 


3.3 


4.7 


5 


279.7 


HI bubble 


IRXS J170849.0-400910 


11.00170(4) 


1.9 


4.6 


8 


742.5 




IE 1841-045 


11.77505(5) 


4.1 


7.0 


7 


529.7 


G27.4+0.0 


AX J1845-0258 


6.9712(1) 






8 


646.1 


G29.6+0.1 



Table 1 . For ea ch of the four AXP s of our sampl e columns 2 and 3 report t he most recent period and perio d derivative from the X-ra y 
timing ('1E1048: Friengo et aljboOSl IRXS J1708: iRea et al.ll2005l IE 1841: iGotthelf et aljl20ol AXJ 1845: iGotthelf fc Vasish3ll998^ . 

column 4 is the surface magnetic field derived from the dipole for mula B = 3.2 X 10i^\ / PP, columns 5 and 6 report the distance 
estimated and the relative dispersion measure derived adopti ng the ICordes fc LaziJ i2002l) model for the distribution of free electrons 
in the ISM (the values obtained usingJ^jfo^^^Cordef ' 199j) model are sim ilar) and last column shows the association with supernova 
remnant or with hydrogen bubble lGaensle^e^aLll20M ; Gaensler et al.l2005^ . Numbers in parentheses are the errors on the last quoted 
digit. 



plored all the DMs up to DMmax indicated in Table 121 For 
the subsequent steps of the analysis (see later), we have on 
the contrary decided to adopt less conservative values for 
the maximum DM explored. 

A typical output plot resulting from the pdm code is 
shown in Figure|5|where the signal of the high-magnetic field 
pulsar PSR J1847-0130 has been searched around its nom- 
inal period and DM in the ranges indicated on the axes of 
the top grey-scale showing the S/N (the darkest the highest) 
as a function of the trial values of P and DM. The central 
panel on the left shows the intensity of the pulsed signal 
as a function of the pulse phase in three sub-integrations 
three minutes long each. The central panel on the right 
shows the folded time series in eight frequency sub-bands. 
The bottom panel is the integrated pulse profile obtained 
summing the sub-integrations and sub-bands using the ap- 
parent (topocentric) P and DM which give the maximum 
S/N. These values, as well as the value of the period as seen 
in the solar system barycenter (BC) and the best S/N are 
reported on the plot. 

For each of our targets we have split the total range of 
the searched parameters in ~ 40 sub-ranges, producing and 
visually inspecting a plot for each of them. A pdm plot has 
been considered interesting - i.e. representing a promising 
radio pulsar candidate - if it showed a high enough S/N 
(above a threshold of 6), a well defined peak in the bottom 
panel and clear linear trends in the grey-scales of sub-bands 
and sub-integrations. 
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Figure 2. Output plot of the pdm programme for the high mag- 
netic field pulsar J1847-0130. From the top, the panels show a 
grey-scale of the S /N as a function of the trial values of P and DM 
(top), the intensity of the pulsed signal as a function of the pulse 
phase in three time sub-integrations (left) and in eight frequency 
bands (right) and the integrated pulse profile (bottom) obtained 
with the P and DM values maximising the S/N. 



With the aim to search for signals from the AXPs with 
periods different from the ones detected in X-rays and in 
the hope to discover new radio pulsars, the data for all 13 
beams have been analysed searching for pulsations over a 
wide range of periods. To do so the pmsearch code (see e.g. 
[Manchester et al. 2001) has been used. This code is based on 
the analysis of the Fourier spectra of the de-dispersed time 
series, obtained, for each different DM, by summing with 
the appropriate delay the 96 frequency channels in which 
the t otal band is splitted. The "Tree" algorithm llTavlorl 
Il974^ used for the de-dispersion works in such a way that 
the maximum DM explored in this case is 10495 pc cm-'', 
corresponding to a broadening of the pulse in a single chan- 
nel equal to the sampling time doubled five times. Pmsearch 
includes also a time domain search algorithm (Fast Folding 



Algorithm; iFaulkner et al .1200411 sensitive to periods greater 
than ~ 3 s, where the frequency resolution of the Fourier 
spectra is very poor and the FFT is consequently less effi- 
cient. 

Pmsearch is one of the blind search codes used for the 
major pulsar surveys performed using the Parkes radio tele- 
scope an d led to the discovery of more than 700 new pulsars 
(see e.g. IFaulkner et alll2004l and references therein). 

Finally a code sensitive to strong single dispersed pulses 
JCordes fc McLauehlinl l20o3) has been applied to the 13 
beams of each pointing with the purpose to discover sig- 
nals similar to those seen in the class of the R otating RAdio 
Transients (RRATs. iMcLauehlin et alj|200d) . The aim was 
to test if AXPs have RRAT-like emission and to search sin- 
gle pulses from nearby sources. 
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3 RESULTS 

3.1 Upper limits on the radio pulsation from the 
observed sample of AXPs 

No radio pulsation has been found in the four AXPs ob- 
served. An upper limit for pulsed radio flux can be estimated 
by using equation0(e.g. Manchester et al. 1996), represent- 
ing the minimum detectable flux density from a pulsar of 
period P. 



sky 



VMHz V 



We 



We 



mJy 



(1) 



Here is the minimum S/N considered (in this case 6.0), 
Tsys and Tsky the system noise temperature and the sky 
temperature in K respectively, G the gain of the radio tele- 
scope in K/Jy, At is the integration time in seconds, A'^p the 
number of polarisations and /S.vmhz the total bandwidth in 
MHz. e is a factor ~ 1.4 accounting for sensitivity reduction 
due to digitisation and other losses. We is the effective width 
of the pulse: 



We 



y/W^ + St^ + 5tl,, + 5tl^, 



(2) 



Its value depends on the intrinsic pulse width W, on the 
sampling time 5t and on the broadening of the pulse intro- 
duced both by the dispersion of the signal in each 3 MHz 
channel (StoM) and by the scattering induced by inhomo- 
geneities in the ISM (Stscatt)- 

Prior to digitisation, the signals are high-pass flltered 
by hardware in order to eliminate slowly varying baseline 
levels. This introduces, depending on the period of the pul- 
sar, a further degrading factor of the sensitivity o f the data 
acquisition system. Following [Manchester et all 1)200 J) we 
estimate a degradation of ~ 40% for periods around 6 — 7 s 
and of ~ 50% for periods around 10 — 12 s. 

Adopting the appropriate numbers for the central beam 
of the multibeam receiver {G — 0.735 Jy/K, Tsys = 23 K, 
Np — 2), an intrinsic pulse width of 5% of the pulse period 
and the values reported in Table|5|we obtain, depending on 
the source and considering the aforementioned degrading 
factors, upper limits in the range ^ 0.06 — 0.1 mJy. 

For 1E1048. 1-5937 the flux density limit would have 
been ~ 0.04 mJy but the observations have been mistakenly 
performed with an offset with respect to the true position 
of the source of 8 arcminutes in right ascension . Assuming a 
Gaus sian beam shape of width 14.4 arcminutes (iHobbs et alJ 
|20o3), a further degrading factor ~ 2.5 caused by the non- 
uniform beam response must hence be added for estimating 
the flux upper limit for this object. 

Folding the time series for IE 1048.1-5937 with a trial 
perio d of 6.16 s and a DM of 334 pc cm~^, the pulsar J1058- 
5957 iKramer et alJ l2003ll has been detected through its 
tenth sub-harmonic, confirming the reliability of the detec- 
tion algorithm used. 



3.2 Discovery of two new radio pulsars 

Analysis of the data collected in all the 13 beams of the 
Parkes multibeam receiver used during this work led to the 
serendipitous discovery of two new radio pulsars, both de- 
tected in the same beam. Since their precise celestial coordi- 
nates are yet to be determined, their provisional names are 



PSRJ1712-3943-1 and PSR J1712-3943-2. The first pulsar 
has a spin period of 0.78 s and a dispersion measure of 525 
pc cm"'', while the second source has a period of 92.5 ms 
and a DM of 713 pc cm"'' (see Table |^. Given the large 
difference in the DMs, we can exclude a physical associa- 
tion between the two pulsars. An immediate confirmation of 
the detection of the new pulsars has been possible thanks to 
the duplication of the pointings performed on three of our 
four targets. The average pulse profiles of the new pulsars, 
obtained summing the data relative to the detection and 
confirmation observations, are shown in Figure |21 

We note that at a distance of 12.2 arcmin from the 
centre of the beam where the new pulsars have been found 
there is the central compact object (CCO) I WGA J1713.4- 
3949 in the supernova remnant G347.3-0.5 jLazendic et alJ 
l2003h . This source could be associated with the 92.5 ms 
pulsar PSR J1712-3943-2: the distance of the CCO from the 
nominal detection position of the radio pulsar is, in fact, 
compatible with the width of the beam of the telescope and 
the DM derived distance of the pulsar D = 7.7 kpc (using 
the lCordes fc Lazij200l model, for which a typical 20% un- 
certainty is expected) is in agree ment with the dist ance of 
the supernova remnant (6 ± 1 kpci lSlane et al.lll999^ . If this 
pulsar is indeed the same object as IWGA J1713. 4-3949, 
we can give a rough estimate of the spin period first deriva- 
tive P that results ~ 10~^* by imposing th at the observed 
X-ray luminosity Lx ~ 6 x lO^^ergs"^ jLazendic et alJ 
l2003h is, as for the bulk of isolated millisecond pulsars de- 
tec ted in X-ray and radio, . 1% of the spin dow n energy 
E (iBecker fc TriimpeJ IT997I: iPossenti et ai]l2002^ . Unfor- 
tunately the period first derivative, calculated using our 
two detections of PSR J1712-3943-2 and assuming that its 
true position is that of the CCO, is totally unconstrained 
{Pobs ~ [1.9±2.7] X 10~^) and cannot be compared with the 
value estimated above. The observations where PSR J1712- 
3943-2 has been detected are infact only few hours apart and 
the periods measured in the two datasets are identical, at l-cr 
level. Using the P derived from the X-ray luminosity we can 
also calculate the putative age of the pulsar r = P/2P ~ 100 
kyr, only a factor of ~ 2 larg er than the age of G347.3-0.5 
( ;S 40 kvr: ISlane et alJfl99g^ ■ a good agreement given the 
crudeness of the estimates. On the other hand, considera- 
tions on the fiux density of PSR J1712-3943-2 weaken the 
probability of an association with IWGA J1713. 4-3949: us- 
ing equatiouQand considering the 12.2 arcmin offset, infact, 
we obtain that PSR J1712-3943-2 would have a flux at 1.4 
GHz 5*1400 ~ 2 mJy at the CCO position. Targetted ra- 
dio observations of the supernova remnant hosting IWGA 
J1713. 4-3949, on the contrary, gave negativ e results provid- 
ing a very tight upper limit of 0.06 mJy iCrawford et alJ 
l2002h . Since radio frequency interferences may have ham- 
pered the detection of the 92.5 ms pulsar in Crawford's et 
al. (2002) single pointing, further observation are neverthe- 
less necessary in order to establish if an association with the 
CCO in G347.3-0.5 (or other close-by high energy sources) 
exists and to determine whether the millisecond pulsar is 
young or recycled. In this second case, given the relatively 
long spin period, PSR J 1712-3943-2 could be a mildly recy- 
cled pulsar. 

The pmsearch code also detected 18 previously known 
pulsars present in the pointed areas. A list of the pulsars 
found during this search (new and known), with detection 
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Name 


P range 
(s) 


(pc cm^'^) 


(pc cm~^) 


r. gal 
"'-scatt 

(ms) 


5t 
(ms) 


(K) 


mjy 


IE 1048.1-5937 


6.0438 - 7.2524 


62754 


682.4 


0.6 


1.2 


5.4 


0.11 


IRXS J170849. 0-400910 


9.6247 - 12.3734 


107069 


1695.5 


127.6 


1.0 


13.1 


0.06 


IE 1841-045 


10.3026 - 13.2451 


114611 


1581.5 


118.9 


1.0 


11.8 


0.06 


AXJ1845-0258 


6.0999 - 7.8426 


67861 


1494.4 


88.8 


1.0 


12.9 


0.06 



Table 2. For each observed AXP we report here the period range and maximum DM explored in the pdmsearch (columns 2 and 3) and 
the values used to calculate (according to equation the upper limit on the flux densities reported in column 8; column 4 reports the 
maxim um DM compatible with the pul sar being in the Galaxy (chosen as the maximum between the values given by the Taylor & CordesI 
il993ft and the lCordes fc Lazid i2002l) models), column 5 the expected broadening induced by the scattering in the ISM at 1.4 GHz, 
column 6 the sampling time and column 7 the sky temperature at 1.4 GHz in the direction of the source. 



Pulsar 


P 


DM 


RAJ 


DECJ 


S/N 






ms 


pc/cm'' 


h m s 


o / // 




J1712- 


-3943-1 


778.149(1) 


525(3) 


17:12:35 


-39:43:14 


16.6 


J1712- 


-3943-2 


92.5277(1) 


713(3) 


17:12:35 


-39:43:14 


12.1 


J1054- 


-5943 


346.9090(4) 


332(2) 


10:54:16 


-59:53:03 


25.8 


J1055- 


-6022 


947.551(3) 


693(7) 


10:56:11 


-60:18:24 


12.8 


J1058- 


-5957 


616.271(1) 


334(4) 


10:58:08 


-59:53:16 


104.4 


J1103- 


-6025 


396.5869(5) 


278(2) 


11:04:03 


-60:18:09 


26.2 


J1705- 


-3936 


854.482(2) 


598(5) 


17:04:58 


-39:43:14 


13.4 


J1705- 


-3950 


318.9408(3) 


202(2) 


17:04:58 


-39:43:14 


47.6 


J1707- 


-4053 


581.017(2) 


348(5) 


17:07:30 


-40:34:00 


10.3 


J1713- 


-3949 


392.4514(5) 


342(2) 


17:12:36 


-39:43:14 


14.4 


J1839- 


-0459 


585.319(1) 


235(3) 


18:39:22 


-04:56:11 


27.8 


J1841- 


-0425 


186.1478(1) 


325(2) 


18:40:21 


-04:31:00 


29.4 


J1841- 


-0524 


455.7312(7) 


283(3) 


18:42:17 


-05:21:23 


12.8 


J1842- 


-0359 


1839.95(1) 


187(9) 


18:41:19 


-04:05:26 


9.7 


J1843- 


-0459 


754.963(2) 


434(4) 


18:43:16 


-04:56:11 


142.6 


J1844- 


-0244 


507.725(1) 


421(3) 


18:44:53 


-02:56:42 


17.1 


J1844- 


-0256 


272.963(1) 


822(4) 


18:44:53 


-02:56:42 


21.0 


J1844- 


-0433 


911.028(3) 


128(6) 


18:44:15 


-04:30:47 


160.0 


J1844- 


-0538 


255.7013(4) 


413(2) 


18:44:16 


-05:21:32 


9.7 


J1845- 


-0316 


207.6358(1) 


503(2) 


18:45:51 


-03:21:53 


31.1 



Table 3. New (first two rows) and redetected pulsars. From left to right pulsar name (for the new pulsars is only a provisional name 
given the position uncertainty), detection period and dispersion measure, right ascension and declination (in coordinates J2000) of the 
centre of the beam containing the pulsar and detection signal-to-noise ratio. The number in parentheses are the 2-sigma errors on the 
last quoted digit. 



period, DM, beam position and S /N (calculated in the time 
domain) is given in Table 

None of the pulsed signals found during this blind search 
could arise from any of the AXPs that we observed, because 
of the discrepancy in position. Figure 2] shows, for four dif- 
ferent values of the DM, the upper limits on the flux density 
at 1.4 GHz reached for the four pointed regions of the sky 
as a function of the spin period. These limits have been es- 
timated using equation with no- = 8 and multiplying the 
limit obtained by a factor ~ 2 taking into account los ses due 
to the search technique (see [Manchester et alJl2'001^ ■ Note 
that the Ua threshold adopted here is higher with respect 
to the one used in the previous step of the analysis: this is 
because, given the huge amount of radio frequency interfer- 
ences (RFI) with which we have to deal when doing a blind 
periodicity search (especially at low S/N), we cannot safely 
state as a good candidate a suspect with S/N less than, usu- 
ally, 8. The RFI problem is not so important when folding 



the time series at a specific value of the period, hence in 
that case we can be less conservative using a lower thresh- 
old for the candidate selection and, consequently, for the rio- 
to adopt in the Smin calculation. 

In order to verify the reliability of the 
quoted flux density limits, we can compare the 
catalogued fluxes of the 18 redetected pulsars 
(obtained from the ATNF Pulsar Catalogue: 
http : //www. atnf . csiro . au/research/pulsar/psrcat/) 
with those estimated putting in equation Q their detection 
periods, DMs and S/Ns (or mean S/N, for pulsars with mul- 
tiple detections) and considering the offset of the detection 
beam with respect to the true position of the pulsars. The 
ratio between the estimated and the catalogued values of 
the flux densities are shown, as a function of the rotational 
period, in flgure |3 Note that the S /N used to estimate 
the flux density, may sensibly vary from one observation to 
another (e.g. because of the different response of different 
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Pulse Phase 



Figure 3. Mean 1374-MHz pulse profiles of the newly discovered 
pulsars PSRJ1712-3943-1 (left) and PSR J1712-3943-2 (right) 
obtained by adding the data of the discovery and confirmation 
observations. The maximum of each profile is placed at phase 0.3 
(the range on the x-axis of each panel being to 1). For each 
profile, the pulsar name, period and DM are given. The small 
horizontal bar drawn under the DM indicates the effective resolu- 
tion of the profile, calculated by adding the bin size, the sampling 
time and the effects of interstellar dispersion in quadrature. 
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Figure 4. Flux density limits for the four observed AXPs, and 
surrounding beams, as a function of the spin period and for dif- 
ferent values of the DM (from bottom to top: 0, 500, 1000 and 
5000). Panels A, B, C and D refer respectively to the pointings 
centered on IE 1048.1-5937 , AX J1845-0258 , IE 1841-045 and 
IRXS J170849.0-400910 . For panel A the solid lines refer to the 
flux limit for the pulsation from the AXP (taking into account the 
8' offset of the central beam with respect to the source position) 
while the dashed line is the flux limit for the blind search in all 
the 13 beams. 



beams in which the pulsar is detected, or because of RFI), 
leading to a large uncertainty in the estimate of the flux 
density. Averaging the S/Ns from different observations 
would lead to a better match between the estimated and 
catalogued fluxes. 

Finally no new RRAT-like signal has been found in 
the single pulse search over the 52 observed regions, while 



4 - 




0.5 1 1.5 

P (s) 



Figure 5. Ratio between the estimated (through equation [Tl 
and the catalogued value of the flux density of the 18 redetected 
pulsars. The errorbars are drawn assuming, besides the errors 
on the measured values, a 20% error on the estimated fluxes. The 
mean value of the ratio calculated over the 18 sources is 1.15±0.3 
(2-(T error). 



two previously known pulsars (J1705-3950 and J1844-0433) 
have been detected also through their single pulses. 



4 DISCUSSION AND CONCLUSIONS 

No radio pulsation has been found from the four southern 
AXPs observed, down to a limit of Smin ~ 0.1 mjy (Table 

Comparing the upper limits on the luminosity at 1400 
MHz (defined as I/1400 = Smin x D^) for our targets with 
the luminosity of the observed radio pulsar population in the 
Galactic field (Figure |SJ, we note that the limits reached by 
our search are a factor of six lower than the median of the 
population. The values obtained for our sample are in fact 
in the range 2.8 - 3.8 mJy kpc^ (shaded region in Fig. El 
while the median of the distribution in luminosity of the 
observed radio pulsars in our Galaxy is 18.0 mJy kpc^. If 
the luminosity function of AXPs refiects that of observed 
canonical pulsars we can calculate that the probability for 
all our four targets to have luminosities below our limits is 
only 5 X 10"*. 

Similar results are obtained comparing the upper limits 
on Li4oo for our sample with the luminosity distribution of 
only young pulsars (with characteristic ages < lO'* yr) or 
long period pulsars (with P > 3 s). 

A comparison with the few high magnetic field radio 
pulsars (having magnetic field strengths above the quan- 
tum critical field line; Figure is not statistically signifi- 
cant, having only four objects in this class. However we note 
that all their lu minosities, as well as the radio luminosity of 
XTEJ1810-197 iCamilo et alJl2006D . are greater than the 
limits obtained for the AXPs. 
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Figure 6. Distribution in luminosity at 1400 MHz of the observed 
pulsar population of the gal actic field (solid line) and of the intrin- 
sic population (dashed line: lLorimer et al.ll200g) . The shaded re- 
gion indicates the luminosity limits reached in the present search. 



197 radio emission is still on more than 3 years after the X- 
ray outburst. Interestingly the sole other possible transient 
AXP is the candidate AX J1845-0258 , one of our targets. 
Our radio observations of this source were performed more 
than six years after its possible X-ray outburst occurred in 
1993, hence unfortunately nothing can be safely concluded 
from our upper limits, in favor or against the possible radio 
and X-ray correlation during the outbursts of this source. 
However, assuming that AX J1845-0258 experienced, after 
the X-ray outburst, a phase of radio emission similar to that 
of XTE J1810-197 , our null detection implies that the fad- 
ing of the radio emission has a time scale of the order of 
few years: in particular, if AX J1845-0258 at the onset of its 
putative radio emission phase had a similar luminosity as 
XTE J1810-197 , this would imply a decrease in I/1400 of a 
factor of ~ 20 over six years. 
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Note however that, if we compare our results with the 
intrinsic lumino sity distributio n of the radio pulsars (dashed 
line in figure IHl iLorimer et al.ll20f)fil) . we obtain a probabil- 
ity of ~ 76% that our observations are not deep enough to 
detect a radio pulsed signal from our targets. 

In considering the causes of the non detection of a ra- 
dio pulsed signal from our four targets, besides the lumi- 
nosity bias, we must take into account the possibility that, 
although the X-ray beam is pointing toward us, the radio 
beam, usually narrower, is not. Assuming a pulse duty cycle 
of ^ 5% (hence a beam semi-aperture ^9°), typical of long 
period pulsars and similar t o that of the radio pulsed sig- 
nal detected bv lCamilo et a l. (2006) in the transient AXP 
XTEJ1810 -197(^ 4% at 1.4 G Hz), we can calculate, fol- 
lowing e.g. iBurgav et alJ ll2003ll . that the probability that 
such a narrow radio beam misses the earth is ^ 77%. The 
composite probability that the beams of all four AXPs are 
not pointing toward us is hence ^ 34%. 

The non detection of RRAT-like bursts from any of 
these AXPs, despite our long exposures, seems to weaken the 
hypothesis that RRAT bursts might be related to the short 
bursts observed from the magnetars leaving us with other 
plausible conjectures of a relation with other classes of neu- 
tron stars such as middle aged radio pulsars displaying giant 
pulses jRevnolds ct al. 2006; Wcltcvrc dc ct al. 2006) o r with 
X-ray Dim Isolated Neutron Stars dPopov et alJl200^ . 

The only case of a detection of radio pulsations from 
an AXP concern s the only confirm ed transient magnetar 
XTEJ1810-197 JCamilo et al.ll2006l) . Radio emission from 
this source is strongly relate d with the occurrenc e of an out- 
burst of its X-ray e mission llHalpern et alJl2005^ . as well as 
an IR enhancement jRea et all2004^ . Furthermore, whereas 
the X-ray fiux is de caying exponentially with timescale of a 
few hundreds days jGotthelf fc Halpernll2005l) . XTE J1810- 
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